function [g2_v, T_order, T] = dynamic_g2(y, x, params, steady_state, T_order, T)
if nargin < 6
    T_order = -1;
    T = NaN(69, 1);
end
[T_order, T] = M30hr.sparse.dynamic_g2_tt(y, x, params, steady_state, T_order, T);
g2_v = NaN(35, 1);
g2_v(1)=T(17)*T(22)-T(18)*T(23);
g2_v(2)=T(4)*T(39)-T(6)*T(40);
g2_v(3)=T(17)*T(36)-T(18)*T(6)*T(35);
g2_v(4)=T(1)*T(43)-T(5)*T(44);
g2_v(5)=T(22)*T(31)-T(23)*T(32);
g2_v(6)=T(1)*T(47)-T(5)*((-y(20))*T(35)*T(48)-T(48));
g2_v(7)=T(4)*(params(4)-1)*params(8)*T(49)-(params(2)*getPowerDeriv(y(15),params(7),2)+T(6)*(params(4)-1)*T(49));
g2_v(8)=T(31)*T(36)-T(32)*T(6)*T(35);
g2_v(9)=T(1)*T(52)-T(5)*T(35)*T(6)*T(35);
g2_v(10)=T(19)*T(26)-params(4)*T(19)*T(27);
g2_v(11)=T(9)*T(53)-T(12)*params(4)*T(53);
g2_v(12)=T(19)*T(38)-params(4)*T(19)*params(1)*T(11)*T(37);
g2_v(13)=T(44);
g2_v(14)=(-y(20))*T(35)*T(48)-T(48);
g2_v(15)=T(10)*T(56)-(1-params(3)+params(4)*T(10))*T(57);
g2_v(16)=T(26)*T(30)-T(27)*params(4)*T(30);
g2_v(17)=T(10)*T(60)-(1-params(3)+params(4)*T(10))*params(1)*((-y(30))*T(37)*T(61)-T(61));
g2_v(18)=T(9)*T(62)-T(12)*params(4)*T(62);
g2_v(19)=T(30)*T(38)-params(4)*T(30)*params(1)*T(11)*T(37);
g2_v(20)=T(35)*T(6)*T(35);
g2_v(21)=T(10)*T(65)-(1-params(3)+params(4)*T(10))*params(1)*T(37)*T(11)*T(37);
g2_v(22)=(-(T(14)*exp(y(17))*T(66)));
g2_v(23)=(-(exp(y(17))*T(28)*T(33)));
g2_v(24)=(-(T(14)*exp(y(17))*T(28)));
g2_v(25)=(-(T(13)*exp(y(17))*T(67)));
g2_v(26)=(-(T(13)*exp(y(17))*T(33)));
g2_v(27)=(-(T(13)*exp(y(17))*T(14)));
g2_v(28)=(-((-1)/(y(15)*y(15))));
g2_v(29)=(-T(49));
g2_v(30)=(-(T(14)*T(68)));
g2_v(31)=(-(T(29)*T(33)));
g2_v(32)=(-(T(15)*T(67)));
g2_v(33)=(-(T(16)*(1-params(4))*T(66)));
g2_v(34)=(-((1-params(4))*T(28)*T(34)));
g2_v(35)=(-(T(13)*(1-params(4))*T(69)));
end
